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@ 
This study was conducted for the State of New Hampshire under Section 22 
Authority. It was performed by Anderson Nichols Co., Inc. under contract to 
the Hydrologic Engineering Section of the Water Control Branch. The study was 
designed to be exploratory in nature and the findings are reported for use and 
comment by field personnel. Responses to the studv report will be used in 
determining the need and direction of any further studies. 
“The main report discusses regression analyses made of recional low flow 
data. Appendix "A" presents the regional low flow data in graphical form 
and Appendix "B" contains summary computer printouts of the “WATSTOR" low 
flow data. 

Supplemental comments and recommendations are as follows: 

1. The developed regression equations indicate zero or negative low flows 
for watershed at or below elevation €00 feet NGVD, and the author suggested © 
(page 59) that, in such cases, a minimum flow value of 0.0001 csm be 
assumed. It is now concluded that there is no basis for selecting this 
value and it is probably excessively low. Since much of the basic data 
used in the analvses was "bunched" about the 1,000 foot elevation (see 
attachments 1 and 2) there is little rationale for extrapolation below 
this elevation. It is therefore recommended that, for watersheds less 
than 1,000 feet in elevation, a minimum "£" of 1,000 feet be adopted for 


use in the developed equations. 


OCLL*= YG POYUCE 


2. It is emphasized that the watershed elevation "E", used in the analyses 
is not the elevation of the stream at the point of interest but is the 
mean elevation of the contributing watershed (see page 39). | 

3. The term, "usable storage", as used in this report, refers to the approx- 
imate regulated storage capacity of a reservoir. Such maximum storage 
capacity should not be interpreted, in all cases, as necessarily usable, 


© as the term might apply. 
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I. INTRODUCTION 


PURPOSE OF STUDY 

Information on low flows is useful for a wide range of planning, 
design and management applications in water supply, hydroelectric 
power generation, fish and wildlife protection, cooling water and 
effluent dilution, navigation, recreation, and the preservation of 
aesthetic features of lakes and streams. With the rapid growth of 
population in. New Hampshire, the demands for low flow information 
have become imperative. The purpose of this study is to investigate 
the magnitude-frequency relationships of naturally occuring low 
flows on streams in New Hampshire to facilitate decision making in 


the management of water resources in the state. 


SCOPE OF STUDY 


The scope of work for this project included the following tasks: 


1) Conduct a literature search to establish the 
availability of any regional low flow frequency studies that may 


have been performed for the State of New Hampshire; 


2) Identify existing water resource projects affecting 
low flows in the respective river basins; ascertain pertinent 
data on these projects and present operating rules and/or 
procedures; assess the probable effect of regulation from 


existing projects on downstream flow; 


3) Develop a simple procedure for estimating low flow for 


natural (unregulated), ungaged streams in New Hampshire. 


In this report, past studies on low flow rates are reviewed. 
The development of a procedure to estimate low flows in ungaged 
watersheds in New Hampshire, and the statistical testing of the 
procedure are described, and a brief list of conclusions is 
provided. Finally, the steps for application of the procedure are 
outlined. The appendices contain low flow frequency curves and flow 
duration curves for the 30 watersheds examined in the study, and 
selected low flow data from the U.S. Geological Survey's WATSTORE 
file for 65 gages. The set of 65 gages is comprised of all New 


Hampshire gages, as well as those gages in Maine and Vermont that 


were included in the statistical analyses for this study. 


e@ 


Il. LITERATURE REVIEW 


INTRODUCTION 

The State of New Hampshire has clearly recognized the need for 
information on low flow discharge rates of streams and for a 
set nodaieay to estimate the low flows of various durations in 
ungaged watersheds. In many other areas of the nation, there is a 
similar need. Recently, the Task Committee on Low Flow Evaluation, 
Methods and Needs, of the Committee on Surface Water Hydrology of 
the Hydraulics Division of the American Society of Civil Engineers, 
conducted an investigation with the following objectives: to 
determine types of low flow information needed for typical problems; 
to describe various low flow characteristics as well as methods and 
data needed for defining those characteristics and their accuracy; 
and to suggest further analyses and data collection needs. During 
its investigation, the Task Committee found that, (1) methods are 
available for defining low flow characteristics from flow records 
and from base flow measurements; (2) methods for estimating low flow 
characteristics at sites without any flow data are inadequate for 
certain purposes; and (3) the use of standardized methods for 
defining various low flow characteristics such as the annual m-day 
low flow or the minimum flow of record, is not possible yet because 
variations in the amount and reliability of data at different sites 
“require that different methods be used depending on the data 
available. This Committee's study did not include an investigation 


of the effects of basin factors on low flows because it believed 


that those effects are "difficult to describe quantitatively and 


therefore are of little value in estimating low flow 
characteristics" (Task Committee, 1980, p. 717). Despite the 
beliefs of the Task Committee, our literature review revealed that 
several investigators have achieved some degree of success in 
analyzing the relationships among low flow discharge rates and 
watershed and climatic parameters. A limited number of these 
Studies have been conducted in northeastern United States. At the 
outset of the present investigation, these past studies were 
reviewed for information on methodology, variables considered in the 
analyses, and findings. A brief discussion of these studies is 


outlined below. 


RESEARCH METHODOLOGIES 
Flow Duration Curves 

As defined by Searcy (1959, p. 1), the flow duration curve is "a 
cumulative frequency curve that shows the percent of time specified 
discharges were equalled or exceeded during a given period. It 
combines in one curve the Flow characteristics of a stream 
throughout the range of sisehates without regard to the sequence of 
occurrence." Flow duration curves are useful in comparing one 
watershed to another. 

The slope of the lower end of a flow duration curve indicates 
the storage characteristics of a watershed: a flat Slope indicates a 


watershed with a relatively large amount of water in storage and 


thus relatively greater low flow rates, while a steep slope 


indicates a watershed with less water in storage and thus relatively 
lesser low flow rates. Estimates of flow duration curves for sites 
on ungaged streams can be made by establishing a relation between 
that stream and nearby gaged streams through plots of concurrent 
-baseflow discharge measurements in the streams. 

A flow duration curve is therefore a summary of all daily flows 
for the period on which it is based, with no regard to sequence. In 
contrast, low flow frequencies, which are based on annual low flows 
for a specified number of consecutive days, are more related to the 
occurrence of drought (Hely and Olmsted, 1963). Low flow frequency 
characteristics of ungaged streams can be estimated from regional 
relations of discharge to drainage area in regions homogeneous with 
respect to geology, topography and climate. The relations are 
developed from concurrent baseflow measurements on gaged streams 
over time. Identification of seepage runs, which consist of 
measuring discharge along a channel reach during a period of base 
flow, can further refine estimates of flows at ungaged sites along a 
stream (Riggs, 1972). In each of these methods, the influence of 


watershed characteristics on low flows is recognized. 


Watershed Characteristics and Low Flows 

The relationship between the physical characteristics of a 
watershed and stream flow has been examined by investigators using a 
variety of methodologies. The factors being studied determine which 
method of analysis best applies. The principal natural factors, 


which are examined individually and in combinations as influences on 


streamflow, include watershed size, precipitation, surficial @ 
geology, elevation, vegetation, and slope. Recent investigations of 
watersheds in New England (e.g., Dingman, 1978; Tasker, 1972; Comer 

and Zimmerman, 1969; Comer and Dunne, 1968; Thomas, 1966) have found 
elevation and surficial geology to be useful in estimating low flow 
discharge rates. 

Low flow discharges are obviously related to precipitation, as 
it is the water input to the physical system of the watershed. The 
relationship between precipitation and low flow discharge rates, 
however, is an indirect one. The volume of precipitation that 
eventually enters a stream channel during a low flow period is 
mediated through several factors: climate, which determines the 
evapotranspiration rates, as well as the amount and form of the 
precipitation; the geology of a catchment, which determines 
infiltration and groundwater flow rates; the topography of the 
watershed, which affects rates of runoff and infiltration; and land 
use and vegetative cover which influence rates of runoff, 
infiltration, evaporation, and transpiration. 

Various precipitation factors have been investigated as 
indicators of low flow discharge rates. Thomas and Benson (1975). 
studied the relationship between streamflow and several watershed 
characteristics in four regions of the United States. They found 
that three low flow parameters ... are significantly related to mean 
annual precipitation but not to precipitation intensity (24 hr, 2 


yr) in the eastern region. However, they cautioned that the 


relations which they defined for low flow estimates are unreliable 


and can provide only rough estimates of low Flow discharge rates at 
ungaged sites. In contrast, Lull and Sopper (1966) found that mean 
daily discharges of the 90% flow duration were correlated with 
precipitation intensity (24 hr, 2 yr) slightly more than with 
average annual precipitation at the nearest weather Statron- 
However, their regression equation is only able to explain 23% of 
the variance in that low flow parameter. 

Other investigators have examined the role of precipitation as 
an influence on low flow in greater detail. Comer and Dunne (1968) 
developed a precipitation index by using recession constants of 
streams in northeastern Vermont to weight daily rainfall amounts 
before and during the times of minimum streamflows. This 
precipitation index is then multiplied by a factor determined from 
either drainage density (for watersheds underlain by poorly drained 
soils), or stream frequency (for watersheds underlain by 
well-drained soils) to predict minimum streamflows. Correlation of 
the precipitation index with the minimum streamflow volumes yielded 
coefficients of determination (R2) ranging from 0.87 to 0.97. 
| Chang and Boyer (1977) used estimates of the September 10-year 
maximum consecutive rainless days and the September mean 7-day 
10-year maximum daily temperature as climatic parameters in an 
investigation of low flows on 12 Monongahela tributaries in West 
Virginia. These climatic parameters were used in the analysis as 
indices of the evaporative power and drought potential of the low 
flow period. This study also demonstrates the value of using a 


precipitation factor which has a Clearly defined and logical 


relation to low flows; the regression equation developed in this 


investigation is able to explain 99% of the variance in the minimum 
7~day 10-year flow. In general, those investigators who have 
attempted to select variables which are representative of the 
physical processes of runoff, infiltration and evapotranspiration 
have been more successful in using a precipitation parameter to 
estimate low flows. 

In many watersheds in New England, streamflow during the summer 
and early autumn consists primarily of groundwater flow. This fact 
suggests that the geology of a watershed may be an important 
influence on low flow discharge rates. As described earlier, Searcy 
(1959) demonstrated that estimates of flow duration curves at sites 
on ungaged streams can be made by establishing a relation between 
that stream and nearby gaged streams through plots of concurrent 
baseflow discharge measurements on the streams. Tasker (1972) 
developed a method for estimating low-flow characteristics of 
ungaged streams in southeastern Massachusetts from basic geologic 
factors and ground water availability maps. He defines a 
"groundwater factor" as a rough approximation of the average 
transmissivity, in hundreds of gallons per day per foot, of the 
watershed. Through regression analysis, the relationships between 
the 7-day mean low flow at the 2- and 10-year recurrence intervals 
and the two independent variables, drainage area and the groundwater 
factor, are defined. 

Hely and Olmsted (1963) believe that the geologic character of a 


drainage basin is the "terrestrial factor" most influential on > 


streamflow. They demonstrated the close relationship between 
lithology and low flow in 19 watersheds in northeastern United 
States. Average values of Yo9/Q, (where Qog = the discharge 
equalled or exceeded 90 percent of time; and Q. = the average 

annual discharge) were determined. For each watershed, the area 
underlain by each geologic formation was weighted according to these 
average values of Y99/Q,- The correlation coefficient between 

the estimates of low flow based on geology and the actual low flow 
For the sub-basins was 0.9. 

Thomas (1966) found a direct relationship between surficial 
geology and low flow in 23 watersheds in Connecticut. Flow duration 
curves for these watersheds show that the "median annual minimum 
7-day average flow ... is 1.30 cfs per square mile from a drainage 
area underlain exclusively by stratified drift, and only 0.013 cfs 
per square mile from an area underlain exclusively by till, a ratio 
of 100 to 1" (Thomas, 1966, p. B209). Streams in the Susquehanna 
River basin exhibit the same relationship between surficial geology 
and low flow values (Ku, Randall, and McNish, 1975). 

Cross (1949) compared flow duration curves of several watersheds 
in Ohio and attributed the differences in shape to variations in the 
geologic characteristics among the watersheds. He concluded that it 
is valid to predict the geologic characteristics of a watershed from 
an analysis of its flow duration curve, but that the converse is not 
a valid approach. | 

A study by Comer and Zimmerman (1969) of two streams in northern 


Vermont attributed differences in discharge rates (of flow durations 


ranging from 1 to 30 days) to differences in topography and soils. To 
Flow duration curves for the two small (3.2 and 8.4 Square miles) 
watersheds were compared. Since the Climate, geology and land use 

are Similar, the infiltration and drainage characteristics of the 

Soil types in the two watersheds are believed to be the major 

influence on the low flow rates of the streams. 

Ives (1977) conducted an empirical analysis of flow duration 
curves in an attempt to define a method for estimating discharge 
rates in ungaged watersheds in New Hampshire. NO Significant 
correlation was found between two baseflow parameters (Q55/0, 
and Yog/Q,) and watershed size for the 14 watersheds in New 
Hampshire. An attempt to relate low flow discharge rates and 
surficial geology produced unsatisfactory results due to lack of 
detailed soils or geologic data. Regional groupings of watersheds 
reflect variations in mean annual runoff; however, this information 
does not aid in the prediction of flow distribution in an ungaged 
watershed. 

In contrast, Dingman (1978) investigated 53 watersheds in New 
Hampshire and found that estimated mean watershed elevation is 
significantly correlated with the ratio of flow exceeded 95% of the 
time to drainage area Qo5/Ap (ft°?/miz). He explains that 
the importance of elevation as an influence on low flow rates is due 
to lower temperatures and evapotranspiration rates, as well as 
greater and longer-lasting snowpack at higher elevations. Thus, 
these factors hold moisture in the watershed and result in a more 


even streamflow throughout the year. Flow duration curves for 


10 


ungaged streams can be approximated by calculating values for the 
2%, 5%, 30% and 95% durations using equations developed in this 
article, and then plotting these values on log-probablility paper 
and drawing a smooth curve through the points. Calculation of the 
95% confidence limits for these flow duration curves is also 
explained. In this investigation, mean watershed elevation is 
estimated by a regression equation developed from data for 
watersheds in New Hampshire. Dingman states that the use of 
measured mean watershed elevation can Significantly reduce the 
confidence intervals about the estimated values. 

The relationship between climate and elevation was also 
recognized by Chang and Boyer (1977) in their investigation in West 
Virginia. They used regressions of elevation and latitude with | 
temperature and precipitation Parameters at weather stations in the 
area, to estimate values of those climatic parameters in the 
watersheds being studied. They pointed out the high correlation 
between elevation and latitude and the temperature and precipitation 
factors. 

Other investigators have attempted to determine the relative 
influence of several watershed characteristics on streamflow rates. 
Lull and Sopper (1966) examined the role of 14 independent variables 
as influences on seasonal average dischage rates as well as on 
various flow durations in 137 small (less than 100 square miles) 
watersheds in the northeastern United States. They found that the 
most influential basin Characteristics are Precipitation, percent of 


the watershed in forest cover, elevation, latitude, average July 


Td 


maximum temperature, and the percent of watershed area that is r 
swamp. However, none of the factors is highly correlated with the 
90% duration flow. 

Studies which examine the relationship between vegetation and 
low flow discharge rates have yielded contradictory results. In two 
studies, forest cover is found to be inversely related to streamflow 
in watersheds in New England and New York (Schneider and Ayer, 1961; 
Lull and Reinhart, 1967). These studies were conducted in 
experimental watersheds and consisted of monitoring streamflow in 
forested catchments before and after clearcutting the basins. In 
each case, streamflow increased after the trees were cut; the 
increases in discharges were attributed to decreased 
evapotranspiration losses. In contrast, Lull and Sopper (1966) 
found a positive correlation between runoff and the proportion of a 
watershed in forest. The soenane suggest that forested land may 
integrate a number of factors, such as slope, latitude, temperature 
and soils, which produce greater runoff. In sedition: they 
emphasize the likelihood that the yield of streams in these 
watersheds would increase if the trees were cut and the 
evapotranspiration losses decreased. 

Few studies of low flows in northeastern United States have 
identified slope as an important influence on low flow discharge 
rates. Valley slope (change in valley bottom altitude over the same 
valley reach divided by valley length), was found to be significant 
in a regression equation for the 7-day two year low flow, but not in 


the equation for the 7-day ten year low flow (Ku, Randall, and 


12 


McNish, 1975). Lull and Sopper (1966) found a correlation 
coefficient of less than +0.316 between main channel slope and 
average annual and average seasonal (four seasons) runoff. 

Other natural characteristics of watersheds have been examined 
in attempts to determine their relationship to streamflow. Examples 
of these are latitude, aspect, average July maximum temperature, 
percentage of swamp in the basin (Lull and Sopper, 1966); valley 
length, sinuousity, valley slope, length“/area, solar radiation 
factor (Ku, Randall and McNish, 1975) and stream density (Hely and 


Olmsted, 1963). 


CONCLUSIONS 

The relationship between the precipitation that falls ona 
watershed and the discharge of the stream draining is an indirect 
one. Natural characteristics of the catchment influence the runoff, 
infiltration, and evapotranspiration processes which operate in the 
watershed. The influences of natural factors on streamflow have 
been studied in sxpevinentat and natural watersheds. In 
northeastern United States, low flows are found to be positively 
correlated with elevation in New Hampshire and the percentage of the 
watershed underlain by stratified drift in Connecticut and inversely 
related to vegetative cover in both New England and New York. To 
effectively manage the water resources in New Hampshire, an accurate 
meted for predicting low flow discharge rates in ungaged watersheds 
would be desirable. Definition of the relationships between 


watershed characteristics and streamflow is a necessary first step. 
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Ili, REGIONAL ANALYSIS OF LOW FLOW DISCHARGE RATES 
UE LARGE RATES 


DESCRIPTION OF PHYSICAL ENVIRONMENT 
Geomorphology 

New Hampshire can be broken into three general classes of 
landforms: the coastal Plain, the lowlands and foot hills, and the 
mountains. The southeastern corner of the state is mainly coastal 
plain, the center of the state is mainly lowlands/foothills and the 
northern portion of the state is mountainous. There is a gradual 
increase in elevation and surface irregularity proceeding from the 
eaieh to the north. Unconsolidated deposits in the mountainous 
region are predominantly composed of thin layers of till with 
pockets of stratified drift. The areas of stratified drift, which 
consists of sands, gravels, and clays, are found in the valleys of 
the major streams. The lowland/foothills region is characterized by 
somewhat thicker layers of till on the slopes with more Frequent and 
larger areal distributions of stratified drift along the region's 
streams. Deposits of marine clays and sands, interspersed with 
till, which are found in the coastal plain, indicate that the area 
was below sea level at one time. 
Climate 

Two factors in the state geography control major aspects of the 
regional climate: topography and proximity to the ocean. Topography 
appears to play a major, complex role, affecting temperature, 
precipitation, moisture storage, and stream discharges. Moist air 


masses tend to enter the state from the southeast, south, southwest 
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and west. As they pass through the state towards the east, @ 
northeast and north, elevation changes in the ground surface force 

the air mass aloft, causing cooling by decompression, and 

precipitation. The areas of highest mean annual precipitation are 

in the northern portion of the state. 

Topography also affects the mean annual temperature within the 
state. Increases in elevation cause a gradual drop in mean 
temperature (about 3.59 F per 1000 ft. in elevation). Therefore, 
the higher elevations tend to be colder, to have a greater portion 
of their precipitation as snow, and to have snowmelt discharges 
throughout greater portions of the year. 

Proximity to the ocean has a moderating effect on the 
temperature variations throughout the year. Sites close to the 
ocean tend toward milder winters and summers, while inland sites 


tend toward colder winters and warmer summers. 


GENERAL PROCESSES THAT CONTROL WATERSHED LOW FLOWS 

To predict low flow,it is necessary first to understand the 
basics of the operating system that produce low flow. Moisture from 
precipitation flows directly to the stream or surface detention 
areas, such as lakes and ponds, or infiltrates the ground surface. 
The infiltrated moisture is capable of flowing through the soil and 
downslope to the nearest stream, but will generally do so only when 
the Beda moisture storage has exceeded its holding capacity. Since 
soil moisture is withdrawn through the evapotranspiration process, 


soils tend to be drier during warmer seasons of the year. 


id 


Therefore, less water is free for streamflow during periods of 
higher evapotranspiration, and precipitation plays a seasonally 
changing role in stream water supply throughout the year. 

The relationship between precipitaton and streamflow is further 
confounded by the movement of water through the ground. The 
potential volume of water which accumulates in and moves through the 
soil and subsurface layers is directly related to the type of 
geologic material which comprises these layers. Porous materials, 
such as sands and gravels, can collect and store large amounts of 
water and release it freely to streams. Flow of water from these 
porous layers to stream channels augments streamflow during times of 
drought. Watersheds drained by porous aquifers generally have more 
moderate flows during extreme wet and dry periods. In contrast, 
solid bedrock and fine-grained sediments, such as clays, permit 
minimal infiltration and movement of water through the subsurface 
layers. Runoff in these watersheds flows directly over the ground 
surface to the stream channels producing a "flashy" stream. Since 
there is limited potential for groundwater contribution to 
streamflow during times of drought, these streams tend to have 
extreme low, as well as high flows. The watershed geology and its 
inbercerattonenins with other geomorphic characteristics of a 
watershed, such as elevation, slope, and surace Storage, therefore, 


influence the severity of any climatically-induced low flow event. 
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LOW FLOW MODELING 
Selection of Modeling Methodology 

One objective of this study is to provide a procedure that can 
be used to estimate low flows in ungaged watersheds with limited 
data on physical parameters. This implies that the procedure must 
be developed from available data and from an understanding of the 
systems operating within these watersheds. 

Since the study will be applied to watersheds with limited data, 
Simulation models can probably be eliminated as a viable approach. 
Measurements of baseflow and investigation of seepage runs were also 
determined to be beyond the limits of the scope of work for this 
study although several sources had advocated concurrent baseflow 
measurements for specific predictions of low flows in ungaged 
watersheds. The problem can be dealt with through some form of a 
statistical model. Statistical models (regression, factor analysis) 
allow the researcher to quantify the relationships between key 
controlling variables and a dependent variable. In this case, the 
objective is to relate certain etimatic and watershed geomorphologic 
variables to low flow rates. 

Application of statistical models creates some significant 
temporal and spatial issues that must be resolved before the model 
is developed. These models are based on a sample of watersheds with 
known flow, climatic, and geomorphologic data. These data are then 
used to explain the interwatershed variation in flow, or in this 
case, low flow. The independent variables used (climate and 


geomorphology) must reflect both the spatial and temporal changes in 
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conditions that influence low flow. For example, the amount of 
precipitation will be critical in determining peak and low flows for 
any given watershed, but the amount of precipitation will vary from 
watershed to watershed (spatial variation). Also, the distribution 
of precipitation throughout a single year (temporal variation) may 
vary from watershed to watershed. Two watersheds with similar 
annual rainfall may have significantly different dry season 
precipitation, resulting in significantly different low flows. The 
selection of variables should reflect a functional relationship with 
the spatial and temporal variations in low flow controls. 

Once these variables have been selected and quantified, they can 
be related to low flows through the use of linear or non-linear 
regression techniques. These techniques will use interwatershed 
variations in both the independent (climatic and geomorphologic) and 
dependent (low flow) yaniebics to create an equation that expresses 
the relationship between the two. Once the equation has been 
computed and tested, it can be used to estimate low flow events in 


ungaged watersheds. 


Selection of Watersheds 

The statistical analyses performed in this study and the 
resulting regression equations utilized streamflow data obtained 
from U.S. Geological Survey (USGS) stream gaging stations located in 
New Hampshire or in Close geographic proximity to New Hampshire. 
The intended use of these equations is to estimate low flow values 


for "natural" conditions on New Hampshire streams. Therefore, the 
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USGS gages selected for use in the regional analysis eee determined ©. 
to have streamflow data representative of natural low flow | 

conditions. The equations were developed and tested using data 

obtained from 30 USGS gaging stations located in New Hampshire, 

Vermont and Maine. Average and minimum low flow data (in cubic feet 

per second per square mile, or csm) for the l-, 7-, 30-, 90-, 183-, 

and 365-day durations for the period of record for each of the 30 

gages are given in Table l. 

The gages utilized in the regional analysis were selected using 
the following procedure. As stated in the scope of work, gages 
first had to meet two criteria: 1) have a contributing drainage 
area of 10 square miles or greater; and 2) possess 15 years or more 
of continuous records. A list of gages was compiled using recent 
USGS Water Resource Data reports for New Hampshire, Vermont, Maine’ 
and Massachusetts. Several discontinued gages which did not appear 
in the recent USGS Water Resources Data reports were added to the 
list. The final list of stream gages which met the criteria for 
drainage area size and length of continuous record included 89 gages 
located in New Hampshire, Vermont, Maine, and Massachusetts. 
Fifty-four of the gages fall within the boundaries of New Hampshire. 

The 89 gages were further investigated to determine which ones 
provided data representative of natural low flow conditions. The 
initial step of this investigation determined which streamflow 
records were affected by regulation using the information contained 
in recent USGS Water Resources Data reports. Generally, the most 


pertinent information contained in these reports was included in the 
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Station description under the category entitled "Remarks". The data @ 
contained in this category included information pertaining to 
upstream projects that affect natural flows at each gaging station. 

For each gage affected by some form of regulation or diversion, 
a list of the flow-affecting upstream projects was developed. 
Because the recent USGS Water Resources Data reports did not provide 
any data relevant to discontinued gages, information on the 
discontinued gages was obtained using USGS Water Supply Papers which 
were compilations of all past gage data (U.S. Department of the 
Interior, 1954, 1964). These compilations contained station 
description data similar in nature to the data found in the recent 
Water Resources Data reports. Information pertaining to upstream 
projects affecting the discharge data at each discontinued gage was 
obtained from this source. Using information compiled from the USGS 
publications, it was determined that streamflow data at 72 of the 89 
gages were affected by regulation to some degree. 

The next step in the process of selecting gages representative 
of natural low flow conditions involved simultaneous data 
collections from two sources, the New Hampshire Water Resources 
Board (NHWRB) and utility companies, to obtain detailed information 
regarding major water projects, such as recreational reservoirs and 
storage pools for hydroelectric generating facilities, which affect 
low flow data. (see Table 2). | 

The goal of this effort was to determine if sufficient data were 
available to quantify the effects of regulation and enable natural 


low flow conditions to be calculated. The results of these 
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investigations determined that the data were not sufficient to @. 
compute natural low flow values at the gages located downstream of | 
these projects. The 49 gages which came under this classification 

were deleted from consideration in the regional analysis. 

The next stage in the selection process was to investigate all 
Gages which, according to the USGS records, were subject to minor 
regulation by various water resources projects. "Some regulation," 
"some diurnal fluctations," and "slight regulation" are examples of 
the descriptions found in the USGS reports. Often the regulation 
was by unnamed mills or ponds. Twenty-two gages were found which 
fell into this classification. Further investigation of these gages 
was conducted using information found in the files maintained by the 
NHWRB. Because there are a large number of dams on New Hampshire 
streams it was deemed infeasible to review the files of every dam 
located upstream of each gage under consideration; therefore, only 
storage dams with pools large enough to appear on USGS topographic 
quadrangle maps were considered to have a significant effect on 
natural low flow conditions. The assumption that storage pools 
which are too small to be detected on USGS quad sheets have an 
insignificant effect on low flow values, appears to be reasonably 
valid, especially when analyzing low flow for durations of 7 days or 
greater. On the basis of this assumption, the streamflow data at 9 
of the 22 gages were judged to be significantly affected by upstream 
regulation and not representative of natural conditions. Further 
investigation of the saute information contained in the NHWRB 


files also indicated that natural streamflow values could not be 
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calculated given the existing data base. Therefore, these 9 gages 
were not included in the regional analysis. The remaining 13 of the 
22 gages which are subject to regulation by minor water projects 
were deemed to be representative of natural low flow conditions and 
were therefore included in the regional analysis. 

Finally, flow duration curves were plotted for the remaining 31 
gages. A comparision of the general shapes of these flow duration 
curves indicated one watershed (located in Massachusetts) which was 
judged to be Gandonestneative and was thus deleted from 
consideration in the regional analysis. Flow duration curves for 
the 30 watersheds in the regional analysis comprise Appendix A. 

Also eon vacned in Appendix A are low flow frequency curves for the 
7-, 30-, 90-, and 183-day durations for each of the 30 gages. 

Comprising the final list of 30 gages with streamflow data 
determined to be reflective of natural low flow conditions were 17 
gages which are not regulated by any projects and 13 gages which are 
affected by insignificant amounts of regulation. A final map check 
was made by comparing the locations of the 30 gages with the 
locations of all significant municipal water projects (wastewater 
treatment plants and water supply intate facilities) in New 
Hampshire (see Tables 3 and 4). 

Wastewater treatment plants were considered significant if they 
had a permit to discharge 1.0 cfs (0.65 mgd) or more. In order to 
comply with water quality standards, the majority of the wastewater 
treatment plants discharge effluent into large streams which provide 


substantial assimilative capacity. Therefore it was assumed that 
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effluent discharges less than 1.0 cfs (0.65 mgd) were insignificant 
in comparision to the streamflow occurring in the recipient 
streams. Water supply projects that regulated a drainage area of 
1.0 square mile or more were considered to be significant. However, 
the influence of water supply projects on low flow values was 
assumed inconsequential if they regulated less than 5% of the 
contributing drainage area at the next downstream gaging station. 
Information pertaining to these municipal water projects was 
compiled using New Hampshire Water Supply and Pollution Control 
Commission (NHWSPCC) publications (NHWSPCC, various dates). 

The review of all significant municipal water projects (water 
supply intake and gaseowates treatment facilities) revealed that 
streamflow data of the gages used in the regional analysis were not 


affected to a significant degree by these projects. 


Selection of Variables 
Dependent Variables: Low Flow Discharge Rates: 

As stated before, one goal of this study was to develop a 
methodology for estimating natural low flows on streams in ungaged 
watersheds in New Hampshire. Therefore, the dependent data set 
consisted of the 7-, 30-, 90-, and 183-day flows at the 2-, 10-, 
20-, and 50-year return intervals for the period of record of the 30 
gaged watersheds in New Hampshire and bordering areas. Discharge 
data for the 30 watersheds to be used in the study were obtained 
from the U.S. Geological Survey (USGS) WATSTORE file (U.S. 


Department of the Interior, 1980). 
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Independent Variables: Watershed Characteristics: 

The independent variables to be used in this study have been 
broken into two major groups: climatic and geomorphologic. These 
variables must reflect interwatershed variations, and must be in a 
form comparable to the low flow data. 

Climate is critical, as it reflects the moisture which comes 
into the watershed and the allocations of that moisture within the 
watershed. The first variable selected to represent climate was 
annual precipitation deviation from the "normal" precipitation 
(CLAD). This variable was selected to express the "wetness" or 
"dryness" of an individual year. Since potential evapotranspiration 
tends to be relatively constant from year to year, this variable 
(CLAD) should generally reflect the tendency for extreme or modest 
low flow events in any group of watersheds or within the record of a 
Single watershed. 

The second variable selected, estimated mean basin elevation, is 
a geomorphologic variable closely related to climate. The 
topography of the state increases in elevation away from the ocean, 
and the elevation of a watershed expresses aspects of temperature, 
orographic rainfall, and snowpack. As a result of these climatic 
factors, less variable flows and relatively higher low flow. 
discharge rates are expected as elevation increases. However, the 
predominance of bedrock in watersheds at higher elevations tends to 
minimize storage of water within the catchment and result in greater 
flow variability. Elevation, therefore, plays a complex role, as it 


incorporates both climatic and geologic factors. 


ID 


Within the wide range of potential geomorphologic variables, 


three others were selected: main channel slope, watershed area, and 
watershed surface storage. Slopes were thought to reflect general 
aquifer conditions (i.e., the higher the slope, the more rapid the 
discharge and the more variable the discharge). Also, steeper 
Sloped areas tend to have more impervious surface materials such as 
till and bedrock, which speed discharge and create more variable 
flows. | 
Watershed area was selected, as the total volume of potential 
discharge will increase with increasing area. A positive 
relationship between area and low flow discharges is expected. 
Watershed surface storage, as measured by the percentage of 
watershed area that is occupied by lakes and ponds, was also 
selected. Increased evaporation losses with increased area of 
surface storage are expected. In general, however, watershed 
storage increases should be positively associated with low flows, 
Since storage will augment flows during dry seasons. 
The following five variables were selected to statistically 
model low flow events: 
1) CLAD 
2) Elevation: 
3) Main channel slope 
4) Watershed area 
5) Watershed surface storage. 
Other watershed characeristics were considered for use in this 
study, but were not included in the final list of independent 


variables due to their interdependence with the other variables, 
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lack of available data, or their perceived relative unimportance as 
influences on low flow discharge rates. For example, data on 
surficial and subsurface geology, which would indicate aquifer 
characteristics, are not available For most aquifers in the state. 
Therefore, a true geologic indicator variable could not be included 
in the analysis, although mean watershed elevation may be an 
indirect geologic indicator. Land use, such as the amount of forest 
cover or impervious area within a watershed, constitutes another 
category of variables which reflect both the geomorphologic and the 
Climatic processes which impact low flows. However, the amount of 
data was limited and the variation in these factors across the 
watersheds in the study area was expected to be small. Climatic 
Factors which were considered but not selected included mean annual 
peeeioleation: precipitation intensity (24-hr. raninfall, 2 yr. 
probability), and July normal temperature. Other investigators 
(e.g., Lull and Sopper, 1966) have found mean annual precipitation 
to be more closely related to mean annual runoff than to low flow 
discharge rates; therefore it was rejected. Precipitation intensity 
has been included in several investigations of low flows, but the 
physical basis for its correlation with low flow discharge is 
unclear. As a variable, July normal temperature would permit 
comparisons of drought potential between different watersheds, but 
not comparisons of flows within a single watershed over time. 
Geographic factors, such as latitude and orientation of the 


watershed were also considered but their influence on low flows is 


an indirect one, and should be expressed in part by the elevation 
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and precipitation variables. Finally, geometric characteristics 
such as watershed shape and length-width ratios were examined but 
not included in the final list of variables because their influence 


on low flow discharge rates was perceived to be minor. 


Data Collection 

In keeping with the scope of work for this investigation, 
published reports and studies were the source of much of the data. 
Most of the watersheds included in this investigation had been 
studied previously; new data was collected only for those watersheds 
and for variables which information was lacking. 

Average annual low flow rates for a range of durations Cli. oe 0s 
15, 30, 60, 90, 120 and 183 days) for the period of record, as well 
as frequency - discharge values corresponding to these durations 
were obtained for the 30 watersheds in the study from the USGS 
WATSTORE data files. All frequencies were computed using a 
Log-Pearson Type III distribution. These data are contained in 
Appendix B. 

The drainage area in square miles of each watershed was obtained 
From recent USGS Water Resources Data reports and USGS Water Supply 
Papers. 

Data on estimated mean watershed elevation were obtained from 
two sources. A recent article by Dingman (1978) contained this 
information for several watersheds. Using USGS topographic maps for 


the remainder of the watersheds, estimated mean watershed elevation 
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0 ee ee 


was Calculated according to the following formula developed by @ 
Dingman: 
E = € : + 0.324 (E -~E ) (1) 
est(i) min(i) max(i) min(i) 
where 
E = estimated mean elevation in watershed 'i' 
est(i) 
E = lowest elevation in watershed ‘i' 
min(i) (the point on the stream for which 
information is desired; 


E = highest elevation in watershed ‘i’, 
max(i) usually located on the drainage divide. 
Main channel slopes were obtained from articles by Benson (1962) 
and LeBlanc (1973) for several watersheds. For the other 
watersheds, values of main channel slope were calculated according 


to the following formula, using USGS topographic quandrangle maps: 


(E(85) (4) - E (15)(4)) (2) 
where Sj) me 
D(i) 
S(i) = slope for watershed 'i' 
E(85)(4) = the elevation of a point 85% of the 
distance from the watershed mouth to the 
endpoint of the main stream in watershed 'i' 
E(15)(4) = the elevation of a point 15% of the distance 
from the watershed mouth to the endpoint of the 
main stream in watershed '‘'i' | 
Dc; = the length of the main stream in watershed 'i' 
(i) 
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The amount of surface storage area, expressed as a percentage of 


watershed area, was obtained for some watersheds from Benson 


(1962). For the other watersheds, it was calculated by 


planimetering on USGS quadrangle maps (1:24,000 and 1:62,500 scales) 


the areas of lakes and ponds within a watershed. These areas were 


converted to a percentage of watershed area, and then 0.5% was added 


to the value for each watershed to prevent values of zero for this 


variable. The formula for computing this variable is as follows: 


Sti) 
where 
St(i) 
L(i) 
P(i) 


B(i) 


(Lea) + Pciy) 


x 100: 6005. Brgy, -(3) 
Bi) 


surface storage in watershed '‘'i' 


lake storage in watershed 'i' visible at a 
scale of 1:24000 or 1:62500 


pond storage in watershed ‘i' visible at a 
scale of 1:24000 or 1:62500 


watershed drainage area. 


Data for the climatic variable, CLAD, were collected from 


National Oceanic and Atmospheric Administration publications (U.S. 


Dept. of Commerce, 1949-1974, 1964). Twenty-one weather stations 


with 30 or more years of records were selected across the study 


area. Thirty years of annual precipitation deviation from normal 


precipitation data were collected for each weather station. 


Values 


of annual precipitation deviation from the normal precipitation 


(CLAD) were calculated according to the following formula: 
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CLAD(n) = Pin) - P (4) 
where 
CLAD(,) = annual precipitation deviation from normal 
precipitation during year n 
Pin) = annual precipitation for year n 
Pp = climatological normal precipitation based on 


the period 1941-1970 

Twenty (20.0) inches were added to each year's value to 
eliminate negative and/or zero values and allow the numbers to be 
transformed to logarithms for a later analysis. For each weather 
station, recurrence intervals for the CLAD + 20.0 (or CLAD*) values 
were Calculated using the USGS plotting position formula (Recurrence 
Interval = (n+l)/m, where n is the number of years (values) and m is 
the order number) (Riggs, 1972). 

Four maps were constructed, one for each of the four frequencies 
(2, 10, 20, and 50 years), according to the following procedure. 
For each weather station, the value of CLAD* for the desired 
frequency was read off the appropriate graph and plotted on a map. 
The 21 values (one for each weather station) were then contoured. 
Input values for the CLAD* variable for the watersheds in the study 
were interpolated from these contour maps (Figures 1-4). 

The streamflow rates, and the morphologic and climatic data 
comprised the overall data set for the study (see Table 5). The 
selection of specific discharge or climatic values for the various 


statistical analyses is outlined in the following section. 
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Application of the Modeling Methodology 


One stated objective of the project was to predict low flows for 
the following periods: 7 days, 30 days, 90 days, and 183 days. This 
was accomplished by relating the previously discussed five variables 
to the measured low flow data. 

The methodology used was the development of a model using 
regression analysis to predict each of the following: 

a) the 2-, 10-, 20- and 50-year seven-day low flows 

b) the 2-, 10-, 20- and 50-year thirty-day low flows 

Cc) the 2-, 10-, 20- and 50-year ninety-day low flows 

d) the 2-, 10-, 20- and 50-year one hundred and 
eighty-three-day low flows. 

In applying regression, it is good practice to have a sample for 
model derivation, as well as a sample for model testing. The test 
sample helps to insure that the derived model does not incorporate 
major predictive errors. Therefore, the sample of 30 watersheds was 
stratified into two groups. The first group consisted of 19 
watersheds and was used for model derivation, while the second group 
of 11 watersheds was used to test the derived models. Figure 5 
shows the location of the 30 gages as well as the major water 
resources projects within the state. 

For the model derivation group of 19 watersheds, values of the 
7-, 30-, 90-, and 183-day low flows at each of the four return 
intervals being investigated were selected. The value of CLAD* for 
each return interval was read from the appropriate map. The 


associated geomorphologic data were also collected for each 
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watershed. The five climatic and geomorphologic variables were 


assumed to be normally distributed and independent, and the 
relationships between the selected variables wete assumed to be 
linear. All flows were converted to cubic feet per second per 
Square mile (csm) to control for the overwhelming influence of 
watershed area on flow magnitude. ‘the var tab ies were related to the 
low flows through the use of 5 rane eeLseenaice multiple linear 
Tegression program of the SPSS (Statistical Package for the Social 


Sciences, 1970). The resultant equations were: 


O.QOO1E + 0.0086 CLAD* - 0.0288 St + 0.0002S 
+ 0.0001 AREA - 0.1395 (5) 


Q7esm 


QO.O001E + 0.0130 CLAD* - 0.0364 St + 0.00025 
+ 0.0001 AREA - 0.1846 6) 


Q30esm 


0.0002E + 0.0239 CLAD* -0.0539 St + 0.00035 : 
+ 0.00004 AREA - 0.2969 (7) 


Q90csm 


Qig3csm = 0.0002E + 0.0489 CLAD* - 0.0870 St + 0.00085 
; + 0.00003 AREA - 0.5308 (8) 
where Q - discharge in cubic feet per second per square 


mile for the duration of n days 


estimated mean basin elevation, 


©) m 
i 
D> 
O 
* 
I tt 


annual precipitation deviation from normal 
precipitation for the recurrence interval being 
investigated + 20.0 

S = (main channel slope, 

St = basin storage, and 


AREA = drainage area, 
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It is interesting to note that elevation and CLAD* explain a 
most of the variance in the data (see Table 6). These variables 


account for the following percentage explanation of variance of low 


flows: 
Q7esm: CLAD* and E explain 66% of the variation 
Q30csm: CLAD* and E explain 71% of the variation 
QoQcsm: CLAD* and E explain 75% of the variation 
Q183csm: CLAD* and E explain 75% of the variation 


The general increase in variance explanation with the increasing 
period of low flow is probably due to the increased sensitivity of 
the climatic variable to long-term droughts. 

Since the variance explanation of these two variables is 
consistently high and accounts for virtually all of the variance 
explanation, linear regression equations were developed using these 


two variables alone: 


Qr 66m =  0.00014E + 0.00780 CLAD* - 0.20959 (9) 
Qz505m = O.00016E + 0.01209 CLAD* - 0.27042 (10) 
Q90c5m = 0.00023E + 0.02231 CLAD* - 0.42554 (11) 
Qigzcsm = 0.00037E + 0.04496 CLAD* - 0.75804 (12) 


The F-test was used to assess whether the inclusion of 
successive independent variables in the regression significantly 
improved the amount of variation in the dependent variable that is 
explained by multiple regression (Till, 1974, p.140). The results 
of the F-tests (see Table 6) indicate that each of the regressions 


is significant at the 0.05 level of confidence. 
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The predictive ability of the equations was assessed through the 
use of a 't' differences test (Till, 1974). This tests the 
hypothesis that the mean difference between predicted and measured 


low flows is equal to zero, and takes the form 


b= Cd: & 0): 57055) 
where 
i= ‘t' value 
d = the mean difference between predicted and 
measured low flows 
SK = s/{7 
Ss the standard deviation of the d values 
i sample Size of the d values 


All of the computed 't* values fell below the critical tabled 
values at the 0.01 level, allowing acceptance of the hypothesis 
(H,) that the mean difference beewsen predicted and measured low 
Flows is not significantly different from 0. The correlations and 
results of the 't' test suggest that the model should be capable, in 
general, of predicting low flows. 

Before the models were accepted as the final results of the 
study, they were applied to the test sample of ll basins. If the 
correlations and 't' values were similar, then the models could be 
considered capable of predicting flows in ungaged areas. The 


results of the analysis were disappointing: 
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Result @ 


Flow Correlation t value 
Q5 0.77 2.54 accept H, 
Q3, Ou77 5.54 reject H, 
Q99 0.82 2.96 reject H. 
0133 | 0.89 3.45 rejecr i 


Although the correlations were similar, the mean differences 
(between what was predicted and what was measured) were 
Significantly greater than 0 for all but the 7-day low flows. With 
these results, a model of this nature cannot be considered reliable 
for application in ungaged basins. 

For simplicity, one of the critical assumptions made in the 
analysis was that variable interrelationships with flow were linear 
in nature. This may not be the case, so scatter plots of the 
relationships between each variable and low flow were compiled. 
These plots indicated that curvilinear relationships did exist 
between the independent variables and low flows. Since this might 
improve the predictive power of che models, a curvilinear approach 
was then used. 

All of the data were reviewed to insure that no negative or zero 
values existed. The climatic, elevation, and flow data were then 
transformed into base ten logarithms, and subjected to a second 
multiple linear regression analysis. As before, stepwise multiple 
regression was performed, in wHtek the variables are entered into 
the equation in order of decreasing contribution to variance 


explanation. The resultant models were: 
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LOG Ql osm = 1.65 Loge 
Log Qsne6m = Ls D2 Log E 
LOD Yognem = 1.39 Log E 


Log 21 83¢sm 1.16 Log E 


These curves were then transformed 


Qreem = (E1E) (erape!-95 
Rope. (Es) (CLADE 2? ) 
Maen, (Ee) Monee") 
Siaseaas: (E>) eran 


+ 1.05 Log CLAD* - 
+ 1.10 Log CLAD* - 
+ 1.28 Log CLAD* - 


+ 1.33 Log CLAD* - 


back into standard 


/ 10 


/ 10 


/ 10°° 


form: 


(14) 
(15) 
(16) 
(17) 


(18) 


(19) 


(20) 


(21) 


The correlations were slightly lower than those of the linear 


models, but the regressions were Significant at the 0.05 level (see 


Table B).: 


Application of the 't' test to each model indicated that 


the differences were not significantly different from 0 at the 0.01 


level. These models were then applied to the test sample of 11 


watersheds, creating the following 


Flow Correlation t v 
Q- 0.80 
Q56 0.79 
oq 0.83 
Q183 0.91 


results: 

alue Result 
2.97 reject H, 
3.01 reject H 
3.08 Teject H, 
4.31 reject H, 
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Again, the model failed to predict values that exhibit little 


difference with measured values. 


A review of the sources of major differences within the 19 


sample revealed an interesting characteristic. Error increased with 
the size of low flow. Since the regression technique used to derive 


these models is based on logarithms, the curve must pass through the 


basin 


origin of the system of data values. If a second regression is 


applied to equations (18), (19), (20), and (21), the curve 


orientation and intercept values can be adjusted to fit the data set 


more effectively. 


Since we are dealing with a bivariate situation (both CLAD* and E— 


are combined as one variable), simple regression was used. The 


resultant equations were: 


Q7csm =i: 23 [e1-65) (cLapel-05) / 107-48 | - 0.0123 
Q30csm = 1.16 | (e1-52)(cLapsl.10) / 106-97] - 0.0134 
teas = 1.10 | (1-39) (cLap*1.28) / 106-56 | - 0.0115 
Qig3csm = 1.17 | (e-16) (cLap*l-33) / 105-62] - 0.0361 


All of the regressions were significant at the 0.05 level, 
the correlations show a marked improvement over the unadjusted 
curves (see Table 6). In all but one case (Qig3), the 
correlations are equal to or better than those of the original 


equations. 


| 


(22) 


(23) 


(24) 


(25) 


and 


Models (22), (23), (24), and (25) were also subjected to 't' 
testing, and the mean differences were not Significantly different 
from O for all four models. The above models were then applied to 


the test sample of 11 basins, with successful results: 


Flow Correlation t value Result 
Q, 0.80 2.05 accept H, 
30 oC 2.45 accept H, 
Qon 0.83 2.5] accept Ho 
Q183 0.91 © 2.69 accept H, 


These results indicate that the derived models are capable of 
predicting flows in ungaged watersheds, will explain up to 80% of 
the flow variance, and should be able to predict low flows of 


intermediate recurrence intervals (2 years<Flow<50 years). 


DISCUSSION OF RESULTS 

The technique used here has provided a reasonable method for low 
Flow estimation. The test sample of 11 watersheds is scattered all 
over the state and adjacent areas of eastern Vermont, and represents 
the range of conditions over which this model will apply. Used — 
within the limits of the sample (watershed size and geographic 
limits of the sample area), this model should be an effective tool 
for estimating low flow discharges. It should be noted that because 
the final equations do not pass through the origin, the derived 


models may give negative flow estimates for streams in watersheds at 
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low elevations (mean watershed elevation less than 600 feet). In @ 
the case of negative flow estimates, a minimum estimate of .Q@001 csm 
could be used. 

There are several areas for potential model improvement. The 
first area deals with data inputs. Because of the limited number of 
watersheds sampled for flow, and the limited physiographic data 
base, a number of more productive techniques had to be eliminated 
from the study. As the data set improves, this should improve 
predictive potentials in this area. 

A second problem is involved with the statistical procedures 
used to derive the final models. There are potential issues 
involved with low flow data serial correlations that could not be 
considered due to the limited amounts of data. As such, these 
potential issues were assumed to be "not significant". As the data 
set improves, there will be an opportunity to research this further. 

The regression techniques used here are not a standard 
polynomial or harmonic approach, and may violate certain assumptions 
of input data normality. The consistent oredictive power of the 
models exhibited in the analysis of test watersheds is judged to be 
indicative that no significant error is introduced by these 
procedures. 

Finally, while the derived curves do not purport to cencesent an 
optimum fit to the data, "optimizing procedures" were utilized which 
have produced a curve that approaches optimum fit to the data. A 


further detailed statistical analysis of residuals would be required 
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@ to significantly improve curve form. However, it is judged that the 
analysis would not significantly improve the oredictive powers of 


the models. 
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IV. CONCLUSIONS 


Four equations (one each for the 7-, 30-, 90- and 183-day 
durations) have been derived which are capable of estimating 
Flows in small to intermediate sized watersheds (10-600 square 
miles). As noted previously, caution must be used when applying 
these models to low-lying (mean elevation less that 600 ft) 
watersheds. 
Fach model is capable of estimating along a continuum of events, 
suggesting that the selected variables are strongly associated 
with interwatershed fluctuations in low flows. 
These equations are valid within the state of New Hampshire. 
Extreme caution must be used when applying these models outside 
the limits of New Hampshire, as variable interrelationships may 


. change. 
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V. APPLICATION OF THE MODEL 


The equations developed in this study can be used to estimate 
flows of 7-, 30-, 90- and 183-day durations, in intermediate sized 


(10 to 600 square miles) watersheds in New Hampshire. The equations 


are: 
1.65 1.05 7.48 
Q osm 1.23 le ) (CLAD* ) / 10 Je 0.0123 
1.63 1.10 6.97| 
Gsnocm 1.16 Ee ) (CLAD* ) / 1o - 0.0134 
1.39 1.28 6.56 
OS aen 1.10 le ) (CLAD* ) / 10 | - 0.0115 
: [1.16 1.33 5-62] 
Q183csm = 1.17 le ) (CLAD* ) / 10 ~ 0.0361 


With maps of contoured values of the climatic variable, CLAD*, 
which are provided in Section III of this report, flows can be 
estimated for four recurrence intervals: 2, 10, 20, and 50 years. 
The procedure for using these maps and the above models is ae 


follows: 


1) Determine size of the drainage area in square miles by 
measuring from USGS topographic maps. 

2) Select the equation for the flow duration period (7, 30, 
90, 183 days) for which flow is to be estimated. 

3) Select the appropriate CLAD* map for the recurrence 


interval (2, 10, 20, 50 years) which is being considered. 
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4) 


>) 


6a) 


OR 6b) 


7) 


Locate the watershed being investigated on the CLAD* map 


and interpolate a value for the CLAD* variable in the 


equation that represents the mean CLAD* for the watershed. 


Using USGS topographic maps, determine the estimated mean 


elevation (E) of the watershed being investigated using 
Equation (1) in this report (p. 39). 

Enter the values for CLAD* and estimated mean watershed 
elevation (E) in the equation and calculate the estimated 
discharge in cubic feet per second per square mile (csm) 
for the selected flow duration and return interval. 

Select the Low Flow Guide Curve for the flow duration for 
which flow is to be estimated (see Figures 6-9). Enter the 
graph with the estimated mean watershed elevation and read 


across to the curve which represents the value for CLAD* 


for the watershed. Then read down to determine the 


estimated discharge in cubic feet per second per square 
mile (csm). 


Multiply by drainage area to determine discharge. 
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APPENDIX A 
FLOW DURATION CURVES 


LOW FLOW FREQUENCY CURVES 


Station No. 


01052500 
01054300 
01057000 


01064300 
01064500 
01073000 
01074500 


01075000 
01076000 
01076500 
01078000 
01084500 
01086000 
01087000 
01089000 
01091000 


01094000 
01130000 


01133000 


01134500 
01135000 


CONTENTS 


Name 
Diamond River near Wentworth, NH 
Ellis River at South Andover, ME 


Little Androscoggin River near 
South Paris, ME 


Ellis River near Jackson, NH 
Saco River near Conway, NH 
Oyster River near Durham, NH 


E. Branch Pemigewasset River 
near Lincoln, NH 


Pemigewasset River at Woodstock, NH 
Baker River near Rumney, NH 
Pemigewasset River near Plymouth, NH 
Smith River near Bristol, NH 

Beards Brook near Hillsboro, NH 
Warner River near Davisville, NH 
Blackwater River near Webster, NH 
Soucook River near Concord, NH 


S. Branch Piscataquoqg River near 
Goffstown, NH 


Souhegan River at Merrimack, NH 


Upper Ammonoosuc River near 
Groveton, NH 


E. Branch Passumpsic River near 
East Haven, VT 


Moose River near Victory, VT 


Moose River near St. Johnsbury, VT 


Station No. 


01137500 


01138000 
01142500 
01145000 
01153500 
01154000 
01155000 
01156000 
01160000 


Name 


Ammonoosuc River near Bethlehem 
Jct., NH 


Ammonoosuc River near Bath, Nb 

Ayers Brook at Randolph, VT 

Mascoma River at West Canaan, Nh 
Williams River at Brockaway Mills, VT 
Saxtons River at Saxtons River, VT 
Cold River at Drewsville, NH 

West River at Newfane, VI 


S. Branch Ashuelot River near 
Marlborough, NH 
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DISCHARGE DATA 
USGS WATSTORE FILE 


CONTENTS 


Station No. : - Name 
01052500 Diamond River near Wentworth, NH 
01054300 Ellis River at South Andover, ME 
01057000 Little Androscoggin River near 
South Paris, ME 

01064300 Ellis River near Jackson, NH 
01064400 Lucy Brook near North Conway, NH 
01064500 Saco River near Conway, NH 

01064800 Cold Brook at South Tamworth, NH 
01065000 Ossipee River at Effingham Falls, NH 
01072850 Mohawk Brook near Center Stratford, NH 
01073000 Oyster River near Durham, NH 
01073600 Dudley Brook near Exeter, NH 
01074500 FE. Branch Pemigewasset River near 

Lincoln, NH 

01075000 | Pemigewasset River at Woodstock, NH 
01075500 Baker River at Wentworth, NH 
01075800 Stevens Brook near Wentworth, NH 
01076000 Baker River near Rumney, NH 
01076500 Pemigewasset River at Plymouth, NH 
01077000 Squam River at Ashland, NH 

01078000 Smith River near Bristol, NH 
01080500 Lake Winnipesaukee Outlet at Lakeport, NH 
01081000 Winnipesaukee River at Tilton, NH 
01081500 Merrimack River at Franklin Jct., NH 


e 01082000 Contoocook River at Peterborough, NH 


Station No. 


01083000 
01084000 
01084500 
01085000 
01085500 


01085800 
01086000 
01087000 
01089000 
01090800 


01091000 
01091500 
01092000 


01093000 
01094000 
01128500 


01129200 


01129500 
01130000 
01135000 
01133000 


01134500 


Name @ 


Nubanusit Brook near Peterborough, Nb 

N. Branch Contoocook River near Antrim, NH 
Beards Brook near Hillsboro, NH 

Contoocook River near Henniker, Nh 


Contoocook River below Hopkinton Dam 
at W. Hopkinton, NH 


Ww. Branch Warner River near Bradford, NH 
Warner River at Davisville, NH 
Blackwater River near Webster, Ni 
Soucook River near Concord, NH 


Piscataquog River below Everett Dam, 
near E. Weare, NH 


S. Branch Pisactaquog River near Goffstown, NH 
Piscataquog River near Goffstown, NH 


Merrimack River near Goffs Falls 
below Manchester, NH 


Sucker Brook at Auburn, NH 
Souhegan River at Merrimack, NH 


Connecticut River at First Connecticut 
take near Pittsburg, NH 


Connecticut River below Indian Stream 
near Pittsburg, NH 


Connecticut River at North Stratford, NH 


Upper Ammonoosuc River near Groveton, NH 


Connecticut River near Dalton, NH 


F. Branch Passumpsic River near 
East Haven, VT 


Moose River at Victory, VI 


Station No. 


01135000 
01137500 
01138000 
01141800 
01142500 
01144500 
01145000 
01150500 
01152500 
01153500 
01154000 
01154500 
01155000 
01156000 
01157000 
01158000 


01158500 
01158600 


01160000 


01161000 


Name 
Moose River at St. Johnsbury, VT 
Ammonoosuc River at Bethlehem Jct., NH 
Ammonoosuc River near Bath, NH 
Mink Brook near Etna, NH 
Ayers Brook at Randolph, VT 
Connecticut River at West Lebanon, NH 
Mascoma River at West Canaan, NH 
Mascoma River at Mascoma, NH 
Sugar River at West Claremont, NH 
Williams River at Brockways Mills, VT 
Saxtons River at Saxtons River, VT 
Connecticut River at North Walpole, NH 
Cold River at Drewsville, NH 
West River at Newfane, VT 
Ashuelot River near Gilsum, NH 


Ashuelot River below Surrey Mt. Dam, 
near Keene, NH 


Otter Brook near Keene, NH 


Otter Brook below Otter Brook Dam, 
near Keene, NH 


~S. Branch Ashuelot River at Webb, near 


Marlborough, NH 
Ashuelot River at Hinsdale, NH 
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